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Quantum fluctuation of the energy is studied for an ultracold gas of interacting fermions trapped in 
a three-dimensional potential. Periodic-orbit theory is explored, and energy fluctuations are studied 
versus particle number for generic regular and chaotic systems, as well for a system defined by a 
harmonic confinement potential. Temperature effects on the energy fluctuations are investigated. 
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A trapped gas of ultracold fermionic atoms constitutes 
an exciting quantum mechanical many-body system that 
presently attracts much attention, see e.g. [H]. The con- 
finement potential, as well as the atomic two-body inter- 
action, can be experimentally controlled |2|, as can the 
number of atoms in the trap. Due to quantum effects the 
energy is expected to vary, or fluctuate, in a non-smooth 
way as the external parameters are varied. A crucial pa- 
rameter in the determination of the fluctuations is the 
dynamics of the system, that may correspond to classical 
regular or chaotic motion. By using semiclassical meth- 
ods 0, 0] energy fluctuations of the many-body system 
are calculated for generic (quantum) chaotic and regu- 
lar systems. We show that a minimum of fluctuations is 
obtained by making the dynamics chaotic, while a har- 
monic confinement potential generally gives rise to much 
larger fluctuations. 

The atoms experience an external trapping potential, 
Vtrap, and a mutual two-body interaction. For a dilute 
gas the atom-atom interaction is well approximated by a 
zero-range delta potential, and the many-body Hamilto- 
nian can be written as, 
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where a is the s-wave scattering length for elastic atom- 
atom collisions, that can be experimentally tuned by a 
magnetic field through a Feshbach resonance [5], from 
large negative to large positive values. We restrict the 
present study to positive values of a, i.e. to repulsive 
interactions between the atoms. In addition, we as- 
sume an ultracold gas, implying quantum effects being 
important. We study a fully unpolarized system con- 
sisting of N fermions with two spin-states, spin-up (|) 
and spin-down (J.). The total density is then given by, 
n (r) = n 1 " (r) (r) = 2n< (r), where n T and are the 
spin-up and spin-down particle densities. From Eq.JT]) 
the single-particle Hartree-Fock equation is obtained, 
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where we have introduced the interaction parameter g = 
4irh 2 a/m. The effective mean field potential is thus given 

by 

eff = gn ] {v) + V trap (r). (3) 



The classical dynamics in this effective potential defines 
the motion as regular, chaotic or mixed and we expect 
different features of the corresponding quantum spectra 
0. The diluteness condition of the gas is fulfilled when 
na 3 <C 1 where nT 1 / 3 is the mean inter-particle spacing. 

With given confinement potential and two-body inter- 
action Eq. Q can be numerically solved, and the to- 
tal energy calculated for a specified number of atoms, 
see e.g. Ref.Q. In this study we are, however, not in- 
terested in a detailed description, but rather in general 
features of the system, and the role of the underlying 
dynamics. This can be obtained by utilizing semiclas- 
sical methods to calculate the total energy of N con- 
fined atoms, U(N) = U + U, at low temperature. As is 
usual in semiclassics [8], the total energy is divided into 
a smoothly varying part U, and a fluctuating part U, 
and we shall focus on the non-trivial fluctuating energy. 
At low temperatures, ksT -C fi, where fi is the Fermi 
energy, the fluctuating part can be calculated as [3j, 
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where the summation p runs over all classical periodic 
orbits in the effective mean-field potential, and r is their 
repetitions. The amplitude A p>r depends on the stability 
of the orbit, S p (E) is the classical action, r p = dS p /dE 
the period and i/ Ptr the Maslov index. Temperature ef- 
fects are included through the function, 
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where tt — h/ (2it 2 kBT). At zero temperature kt = 1. 

The classical functions in Eq. ([4]) are evaluated at E = 
/i. Since Eq.Q suffers from severe convergence prob- 
lems, we characterize the fluctuating energy, U(N), by 
computing its moments. The first non-trivial moment is 

the variance, er 2 = (^U 2 ^ that gives the typical size of the 
fluctuations. This can be written as [3] 
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where K (r) is the diagonal part of the spectral form fac- 
tor (the Fourier transform of the two-point energy corre- 
lation function). The form factor is by definition system 
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dependent but for long times, i.e. r 3> r m i„ (r m , n is the 
period of the shortest periodic orbit), general statistical 
properties can be derived 0, OH] , namely 



for the regular case, and 
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for the chaotic case. In these expressions th — hp is 
the Heisenberg time and 6 the Heaviside step function. 
The so-called r m ;„-approximation assumes the smooth 
behavior of K (Eqs.© and ©) all the way down to 
t = Tmin\ for t < T m i n there are no periodic orbits and 
the form factor is zero. In this approximation simple 
expressions can be derived for the energy fluctuations 
assuming regular or chaotic dynamics, by inserting the 
corresponding expressions for the form factor into Eq. © . 
At temperature zero this gives [3], 
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for the regular case, where E c — h/T min and b = E c p = 
th I Tmin ("dimensionless conductance"). In the range 
Tmin <C t <C Th , and for not too small numbers of par- 
ticles, it is a good approximation to write K c h = It. 
In this approximation the energy fluctuations for chaotic 
dynamics become, 
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In the studies below we shall, however, generally inte- 
grate Eq.© numerically. The important parameters to 
determine the fluctuations in energy are thus the period 
of the shortest periodic orbit, T m i n , and the mean level 
spacing at the fermi energy, S — l/p = dp/dN. As 
shortest orbit we assume a "diameter orbit", i.e. mo- 
tion along the diameter in the effective potential, and 
get T m in = 4 L max dr/v{r), where v(r) is the classical 
velocity and R ma x corresponding maximal radius. 

At non-zero temperature Eq.© is solved numeri- 
cally using the respective expression for the form factor, 
Eqs.© or ([8]). In the limit of very low temperature an 
analytical expression of the energy fluctuations may be 
derived for the chaotic case (tt >> th » T m i n ) [111 ]. 
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and similarly for the regular case, 
o h 2 r H ( 1 



When the trapping potential is harmonic, V tra p = 
mw 2 r 2 /2 (to is the atom mass; we shall put to=1), 



the rather weak two-body interaction between atoms im- 
plies that the effective potential Eq. l(3j) can be approxi- 
mated by uj e ffr 2 /2 + er 4 /4: [12], implying regular dynam- 
ics. Only few periodic orbits appear, and the expression 
Eq.©, describing a generic regular system, does not ap- 
ply. Instead, fluctuations in energy is well described by 
two periodic orbits, namely the circular and diameter pe- 
riodic orbits [12] • The interference of the two orbits gives 
rise to supershell structure, as noticed in 0|. The semi- 
classical expression of the energy in terms of these two 
periodic orbits reproduces very well the result obtained 
from a microscopic Hartree-Fock calculation. 

The second moment of the energy is found by squar- 
ing Eq.©, including only the two shortest periodic or- 
bits. Only the diagonal terms and those corresponding to 
cosines of the action difference contribute, since the oth- 
ers are eliminated through the averaging procedure due 
to the rapid fluctuations in the arguments. The resulting 
expression for the energy variance becomes, 
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where A d and A c are amplitudes of the diameter and 
circle orbits, respectively [12]. The temperature depen- 
dence simply appears in the pre-factor kt (Eq.©). The 
first two terms, corresponding to squares of the two peri- 
odic orbits, diverge at the bifurcation point e = 0. Con- 
vergence of the expression is, however, restored by the 
cross-term. 

In FigQ] we show the energy fluctuations, 02 po , as a 
function of iV 1 / 3 for some different interaction strengths, 
.9 = (pure harmonic oscillator; H.O.), g — 0.1, 0.2 
and 0.4. The energy is here (and throughout the paper) 
expressed in units of hu. Energy fluctuations are seen to 
decrease with increasing value of g. Supershell structure 
implies local minima close to <J2 P o — 0. For larger values 
of g the supernode appears at smaller particle numbers. 
In the case g — 0.4 also the second supernode is seen in 
FigQ] For this interaction strength we also compare to 
the Hartree-Fock result, that is seen to give quite similar 
result. 

Energy fluctuations in a generic regular system is de- 
scribed by Eq.©. Assuming same scaling with particle 
number as for the H.O. for Ec and p, gives energy fluc- 
tuations as shown in FigUl No supershell structure ap- 
pears, but a reg is found to increase linearly with N 1 / 3 . 
The effect of the interaction strength is indeed quite 
small. Compared to the harmonic confinement potential, 
FigQ] the size of the shell energy fluctuations is notably 
smaller. 

A drastic decrease of energy fluctuations thus appears 
as the number of contributing periodic orbits increases 
from one (family) orbit for the pure oscillator (g=0 case 
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FIG. 1: (Color online) Energy fluctuations, <T2 po , versus parti- 
cle number, N 1 ^ 3 , for the H.O. confinement potential with the 
two-body interaction strength g = 0, 0.1, 0.2, 0.4 (black dot- 
ted, blue dashed, green solid and red dashed-dotted). Corre- 
sponding Hartree-Fock result for <?=0.2 is shown by the thick 
solid black line. 
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FIG. 2: (Color online) Energy fluctuations versus particle 
number, 7V 1//3 , assuming a general regular (left-hand figure) 
or chaotic (right-hand figure) effective potential with the two- 
body interaction strength g = 0.1, 0.2, 0.4 (blue dashed, green 
solid and red dashed-dotted). 



in FigHJ), to two dominant orbits for the harmonic con- 
finement with an interacting gas (g > in FigQ]), and 
finally to the generic regular system (Fig(2k; e.g. a po- 
tential with steep walls ) with several orbits contribut- 
ing. The diminishing of energy fluctuations with increas- 
ing number of periodic orbits occurs due to destructive 
interferences from several orbits. The special shell struc- 
ture with large energy gaps and degeneracies shown by 
the harmonic oscillator quantum spectrum, is less pro- 
nounced for systems described by many periodic orbits. 

Classical dynamics in the effective potential may be- 
come chaotic by properly arranging the confinement po- 
tential. In the purely chaotic case the number of periodic 
orbits increases exponentially with the period time of the 
orbit and the T m i„-approximation works well. The fluc- 
tuations of the energy originating from chaotic dynamics 
can thus be calculated by inserting Eq.® into Eq.©. In 
Figl2b energy fluctuations for chaotic dynamics is shown 
versus particle number for different values of the inter- 



action strength. For large particle numbers, N > 200, 
Eq. iflOj) is a very good approximation of a c h ■ 

In the case of chaotic dynamics, the energy fluctuations 
are indeed very small. As compared to the regular case 
with a harmonic confinement potential, see Fig. Q] the 
fluctuations are about four orders of magnitudes smaller 
in the chaotic case. The fluctuations are also seen to be 
fairly independent of particle number and on the inter- 
action strength for N > 100. The drop in a ch at smaller 
particle numbers originates from higher order terms in 
the expression for K c h, see Eq.®. 

Experimental conditions always implies a non-zero 
temperature, and we shall now investigate the temper- 
ature dependence of the energy fluctuations. We assume 
that the system can be described by the same density 
of states at T > as for T = 0. This implies that the 
fermi energy, the periods of the periodic orbits and the 
mean level spacing are all assumed independent of tem- 
perature. This turns out to be a good approximation for 
temperatures here studied. The change in temperature 
thus appears only in the function kt, Eq.®. 

The fluctuation in energy at non-zero temperature is 
obtained by numerically solving Eq. © with the form fac- 
tor given by Eq. ((TJ) for the regular case, and by Eq.® 
for the chaotic case. In the case of a harmonic trapping 
potential the temperature dependence of the energy fluc- 
tuation is directly given by kt, see Eq. lfTTj) . 
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FIG. 3: (Color online) Energy fluctuations versus particle 
number at different temperatures, assuming a harmonic trap- 
ping potential with interaction g — 0.2. The three curves are 
valid for temperatures fcsT/ftw = 0.0, 0.1, 0.2 (blue dashed, 
green solid and red dashed-dotted). 

In Fig[3] the energy fluctuation is shown versus par- 
ticle number for three different temperatures. For the 
harmonic trapping potential with g — 0.2 the supershell 
structure is clearly preserved but is suppressed with in- 
creasing temperature. 

The temperature dependence of the energy fluctua- 
tions is shown in Fig|4]for the three studied cases, reg- 
ular dynamics in harmonic trap (contribution from two 
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FIG. 4: (Color online) Temperature dependence of energy 
fluctuation relative to zero temperature for g — 0.2. H.O. 
confinement (blue dashed line), generic regular in T m i„ (green 
solid line) and chaotic (red dashed-dotted line). Temperature 
is expressed in units of h\o. 

periodic orbits), generic regular dynamics, and chaotic 
dynamics. The fluctuations are shown relative the T = 
fluctuations. All three cases show a similar temperature 
dependence. At a temperature ksT « O.IHuj the fluctua- 
tions have decreased by 50%, and more or less no energy 
fluctuations appear for fc^T > O.Shuj. 

In summary, effects from regular and chaotic dynamics 
have been studied for a trapped dilute gas of interacting 
fermionic atoms. The fluctuation of the total energy was 
studied using periodic orbit theory, and three different 
dynamical systems were considered, regular dynamics in 
a harmonic trap, general regular dynamics, and chaotic 
dynamics. 

In the case of a harmonic trapping potential a very 
small number of periodic orbits are sufficient to describe 
the fluctuating part of the energy. If the atoms are 
not interacting the Hamiltonian has SU(3) symmetry 
(pure harmonic oscillator), and one (family) orbit ap- 



pears. This gives rise to a monotonically increase of en- 
ergy fluctuations with particle number. A week repulsive 
interaction between the atoms implies the contribution 
of two dominating periodic orbits. The two orbits in- 
terfere and give rise to supershell structure, where the 
detailed behavior depends on the interaction strength, g. 
In general, the fluctuations were found to decrease with 
interaction strength g, and to increase non-monotonically 
with particle number. 

For a general regular system, as for example an effec- 
tive potential with a flat bottom and steep walls, several 
periodic orbits contribute. The fluctuations of the en- 
ergy were then found to be proportional to N 1 / 3 , and 
to be considerably smaller than for the harmonic trap- 
ping potential. The fluctuations are not sensitive to the 
interaction strength. 

If the effective potential is arranged in such a way that 
the dynamics is chaotic, the energy fluctuations are found 
to be considerably diminished as compared to the case of 
regular dynamics. Compared to the case of a harmonic 
trap, the fluctuations for chaotic dynamics were found to 
be almost four orders of magnitude smaller for particle 
numbers N > 10 4 . This means that the energy fluctu- 
ations more or less vanish for the chaotic effective po- 
tential, independent of interaction strength and particle 
number. 

Finally, effects of temperature were described in peri- 
odic orbit theory. As expected, a non-zero temperature 
diminishes the energy fluctuations. Similar dependence 
on temperature was found for the studied different dy- 
namical situations, and the energy fluctuation was found 
to be 50% smaller at temperature fcsT w 0.1 in units 
of hu>, as compared to zero temperature. For a trap fre- 
quency of 10 4 s _1 temperatures of the order of 10 nK or 
smaller are thus required to observe significant fluctua- 
tions and supershell structure. 

We are grateful to Stephanie Reimann for useful com- 
ments. 
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